LATTICE SYSTEM OF INTERACTING SPINS IN THE THERMODYNAMICAL 

LIMIT 



S. SERGEEV 

Abstract. In this paper we investigate some particular spin lattice (a higher dimensional generalization of 
a spin chain) related to Zamolodchikov model, in the limit when both sizes of the lattice tend to infinity. An 
infinite set of bilinear equations, describing a distribution of eigenvalues of infinite set of mutually commuting 
operators, is derived. The distribution for the maximal eigenvalues is obtained explicitly. The way to obtain 
the excitations is discussed. 



Introduction 

Integrability of Zamolodchikov-Baxter three dimensional spin model |2 ^ is based on the existence of 
commutative set of transfer-matrices T{6\, 02, 63), 



(1) 



T(8 1 ,9 2 ,8 3 ),T(8 1 ,9' 2 ,8 3 ) 



0. 



where 6j are Zamolodchikov' dihedral angles, and we understand T as an operator in the vertex formula- 
tion [3] of Zamolodchikov-Bazhanov-Baxter model [I| with two spin states. In contrast to two-dimensional 
integrable models, T is layer-to-layer transfer matrix joj, that means 

• it is associated with a rectangular lattice with the size N x M, therefore matrix T has the dimension 
2 NM x 2 NM 

• two parameters 62 and 63 are varied in eq. Q . 

Matrix elements of i?-matrix of the Zamolodchikov-Baxter model are not positively defined: it is the ob- 
stacle for the decent interpretation of it as the model of statistical mechanics. The quantum mechanical 
interpretation is preferable. 

Relation implies the existence of a set of commutative operators {t TO , n (0i)} such that 

(2) [*m,n(0i), T(B X , 8 2 ,8 Z )] = Vff 2 ,fl3,m,n, 

i.e. the problem of diagonalization of T for any 82,83 and the problem of simultaneous diagonalization of 
{t m .n} are equivalent. In what follows, we mean a determined definition of {t m n \ with < m < M and 
< n < N related to an auxiliary problem of the model. 

It is well known, the Zamolodchikov model and its generalization - Bazhanov-Baxter model @] - are 
related to the generalized chiral Potts model 0. The set of {t m , n (8i)} may be produced by the expansion of 
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transfer matrices of the length = M chain of Lax operators for cyclic representation of U q= -i(sIn) (these N 
and M are exactly the size of the layer). Another scheme to produce the same set {t rn . n (9i)} was proposed 
in |Sj ■ This scheme is the invariant one from the point of view of 2 + 1 dimensional integrability, in 
particular the N «-> M symmetry is evident in this scheme. All t m ^ n are simple polynomials in the algebra 
of observables, and besides they are evidently hermitian (i.e. the model is indeed a model of quantum 
mechanics). The invariant scheme we mention here as the system of interacting spins on a two dimensional 
lattice, or as the two dimensional spin lattice. 

The eigenvalues of the set {i mj „} may be found as a solution of a system of second order equations. In the 
language of auxiliary transfer matrices for generalized chiral Potts, the system of second order equation is 
the complete set of fusion relations for fundamental transfer matrices. The reader may find the investigation 
and discussion of the fusion relations and Bethe Ansatz for Zamolodchikov model for TV = 3 in [51E1E3]- In 
the direct 3D scheme the whole system of fusion relations is encoded into a single spectral equation [7J . 

The physical problem is to find the spectrum of all t m ,n(^i) simultaneously. One way to solve this 
problem is the nested Bethe Ansatz for the fusion algebra, M — > oo and finite N. Suppose for a moment, 
somebody has succeeded in solving the nested Bethe Ansatz equation for arbitrary N (i.e. in finding the 
limiting M — > oo densities of N — 1 Bethe Ansatz 's distributions of zeros) and then sends N — > oo. From 

2 + 1 dimensional point of view such 1 + 1 dimensional result would be related to the limit AT, M — > oo with 

N 

— — > 0, i.e. N «-> M symmetry would be lost - the model in this approach remains the 1 + 1 dimensional 
model with infinite symmetry group. This approach would give a correct answer for a quantity independent 
on N/M. 

Contrary to this, the spectral equation in the direct 2 + 1 scheme is initially N <-> M invariant. In this 
paper the spectral equation is evaluated in the limit 

N 

(3) N,M -> oo, — C 

where ( in the non-singular aspect ratio for the layer. The main result of this paper is the exact distribution 
of the largest eigenvalues (the ground state) t m>n = f(m,n;9i,() in the limit ©. The other result is the 
limiting form of the spectral equation allowing one to describe (at least qualitatively) the gap-less excitations 
of the ground state. 

This paper is organized as follows. In sections 1,2 and 3 we formulate first the system of interacting spins, 
recall its finite N X M - volume spectral equation and make its leading term evaluation. Content of the first 
three sections is a repetition of |7I IH1 1111 1121 113j . Next, in the fourth section, we expose some preliminary 
numerical results for the spectrum of t m ,n and discuss the main idea for the limiting @ procedure. In the 
fifth section re- write the spectral equations in the thermodynamic limit N,M — > oo. In the sixth section 
the qualitative analysis of the thermodynamical spectral equation is given, the distribution of the maximal 
eigenvalues of i m ,„ is obtained and the structure of excitations is discussed. 
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1. Formulation of the spin lattice system 

All the ways to produce the set {t m ,ri}, both via Lax operators for cyclic representation of £/ g= _i(s/Ar) 
and via 31? linear problem 0, finally may be reformulated in the following combinatorial form. 

Consider a square lattice with the size N x M with periodical boundary conditions - exactly the layer of 
||TJ. Each vertex j of the lattice may be labelled by the pair of the indices j — (n, m), n € Z^r, m € %m- A 
local triplet of the Pauli matrices m , cr^.m an d a n m = ' a n m a n m i s assigned to each vertex. 

Consider a set of non-self-intersecting paths on the periodic lattice with the following rules of bypassing a 
vertex and following factors jj associated with each variant of bypassing (note the multiplier k in the third 
variant): 




n n n 

7n,m = ^n,™ 7n,m = °n,m 7n,m = K(J n,m 



An example of such path for 4x4 lattice is drawn below: 




Any path V has a homotopy class c{V) = mA + nB, where A is the cycle left to right and B is the cycle 
from bottom to top. In the other words, m is the horizontal winding number and n is the vertical winding 
number of the path V '. The path at the example above has n = m = 1. 
For fixed winding numbers n and m let 

(4) 4,»w = 

V : c(V)=mA+nB V 

be the sum of the products Y\ Ij °f 7-factors along a path V for all possible paths with the given winding 

V 

numbers. In particular, J^o = 1- The winding numbers of J m , n run < m < M and < n < N . The reader 
should distinguish the periodical discrete coordinates (n,m) £ (Zjv,Zm) of the algebra of observables and 
the winding numbers (m,n) labelling the operators J. It is known 0|H|, operators J m ,n obey the following 
exchange relations: 

/r\ IT,, — (_\nm'+n'm j j 

V u y ° m ,n ° m ,n \ } ,J m ,n ° m ,n • 
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It means, they can be quasi-diagonalized simultaneously: 

(6) J m M = (^) m K) n * m ,n(K) , [**,,«(«)»*»»',»'(«)] = 0. 

Auxiliary matrices er^ and er a , cr^cr^ = —a v a x , belong to the algebra of observables. They correspond to the 
homotopy classes 1 • A + • B and • A + 1 • B. Without loss of generality one may fix 

(7) a* = []<i and <r» = 

n m 

The key meaning of the auxiliary matrices is that they represent one extra degree of freedom of {J m ,n} with 
respect to {t m n }. For any of 2 NM ~ 1 eigenstates of {i mi „} auxiliary a x ,a v are usual 2x2 Pauli matrices. 

The set of J m>n (and {t m n } as well) is the set of "integrals of motion" for Zamolodchikov model [2] 
in its vertex formulation [3] . Namely [7J |SJ I12j , the layer-to-layer transfer matrix of Zamolodchikov model 

Q 

T{0±, 62, 83) commutes with all J m ,n for « = tan — and arbitrary 62, 0%. In this paper we prefer to call t m ,„ 
the moduli since in the classical limit they become the moduli of the classical spectral curve |14| . 

The advantage of the present quantum-mechanical formulation is that if k is real, all J m>n and t m>n are 
self-adjoint since the Pauli matrices are self-adjoint, therefore the model is evidently physical. An eigenstate 
of the model is defined by eigenvalues of all t m ,n - one can label the eigenstates by the corresponding values 
of {t m ,n}- Our aim is to describe all eigenstates. 

2. Finite size spectral equations 

In this section we recall the functional equation for the set of t m ,n- For its rigorous derivation see |1 1| . 
Consider the following generating function: 

M N 

(8) j(x, y ) = J2 E(-)" +m+ " ma;m y" J ™.™ > 

m— n— 

where x and y are generic complex numbers. In the basis of the auxiliary er-matrices J(x, y) is 

(9) J(x,y) = to,o(x,y) - o- x tifl(x,y) - a y t 0yl (x,y) - a z t 1:1 (x,y) 
where 

(10) 

to,ifav) = E(-)" x2m y 2 " +1< 2™,2«+i , t hl (x, y ) = ]r(-r +? v™+y" +i t 2m+1 , 2 „ +1 . 

rn.n m,n 

The reader should not be confused by the notation t a p(x,y) with a, — 0,1 and the set of t m n , < 
m,n < M,N. It is known ^203^], the complete Abelian algebra of t m ,n is generated by the polynomial 
decomposition of 



(11) 



to,o( x -.y) 2 ~ hAz^y) 2 - to,i(x,y) 2 - t 1:1 (x,y) 2 = F(x 2 ,y 2 ), 
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where 

N-l M-l 

(12) F(X N ,fX M ) = ]"( (1 - \e 2 ™ n / N - ^"Kim/M _ K 2^ e 2wi(n/N+m/M)j 



n— m— 
lJV _ _2 ,,M _ „,2. 



is a polynomial of A = x , /r = y : 

M N 

(is) iV, y 2 ) = £]>> 2 V q ^q. 

P=0Q=0 

As it was mentioned in the introduction, equation l|lll) encodes the whole fusion algebra of auxiliary transfer 
matrices for U q= -\ {sis), the reader may find its explanation for e.g. N = 3 in the appendix. 
The right hand side of may be re-written as 

(u) J2 x ' 2P y 2Q T,T,(-r +n+mn+mQ+nP t m .nt2P- m .2 Q ^ ■ 

P,Q m,n 

Equation is the principal solution of the model, in the same way as the Bethe-ansatz is the principal 
solution for the spin chains: the problem of diagonalization of NM 2 x 2 NM matrices t„ l; „ is reduced to 
a system of NM algebraic equations. 



(15) ^^(-) m+n+ " m+mQ+nP i m: „i 2 p- m ,2Q-, i = F P . 



Q 



3. The leading term and relation to Zamolodchikov model 

Suppose, no one term is zero in the product (|12fl . Then F in is exponentially big, and one may 
definitely conclude [T3). 

(16) Each of (t a ,ptx,y)) ~ \F{x\y 2 )\ ~ e™«(W) , 

V /a, 0=0.1 



where x = A^ 2 , y = l^ M ^ 2 , and the integral 



i i r r 2 * , 

1 , 1 / / ,,,,,, . ;^ 2 X((p i(0+0')| 



(17) fl(A lA1 ;« 2 ) = lim — log |F(X, F)| = — w / / log 1 1 - Ae 1 * - /xe'* - K 2 A/xe 



being parameterized by 



sin r 2 , , sm r 3 2 sin rp sm ri 



(18) |A| = — - , |/x| = — — , = 

sm n sm n sm r 2 sm r3 

with r,- bounded by 

(19) r + n + r 2 + r 3 = 7r and < n + r 2 , n + r 3 , r 2 + r 3 < 7r 
has the value ^3] 

3 

(20) q{\,h;k 2 ) = -log2sinri (— log 2 sin , 
where <£>(r) is the polylogarithm pQ: 

(21) *(r) = f] Sin(2mr) 



2 



27rm 

771=1 

This value is closely related to Baxter's result for the bulk free energy of the Zamolodchikov model [Tl ll2l[T^] . 
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Relation l|ltj|) gives the solution of I|1U|) for the mean eigenvalue problem for t a: g{x,y): any eigenvalue of 
t a ^(x,y)^j has this leading behavior. The answer i|20|) is unsatisfactory from the quantum mechanical 
point of view, it corresponds to the asymptotically infinite values of the spectral parameters x — \ N / 2 and 
y = fx M ' 2 and do not clarify the structure of eigenvalues of the whole set of i m ,„. 

4. Preliminary evaluation for finite N, M 

We started the investigation of H15|) with the numerical tests for relatively small N, M (up to N = M = 8) 
and for simple choices of k (k = 0, 1). 

The principal observation for finite N, M is the following. Excluding i m ,„ from (|15fl step-by-step, one 
comes to a final polynomial equation for a single t m n : such polynomial equation is exactly the characteristic 
equation for the operator t rn n . Therefore, the system (|15|l and the problem of direct diagonalization of 
operators t m ,„ are equivalent. In other words, any solution of equations (|15H is indeed an eigenstate. For 
this reason we call (|15H the complete Abelian algebra. 

Note in addition the parity property: if a set {t m ^ n } solves the equation l|15fl. then the set {t rn . n }, 

(22) t2m+a,2n+0 = £a,/3^2m+a,2n+/3 , 

where a,0 = 0, 1 and e a g are four arbitrary signs, solves l|15f) as well. This ambiguity corresponds to the 
ambiguity of definition of the auxiliary <j x ,a v ,a z . 

It is useful to visualize the domain of the indices of {t m n } as the set II of points (m, n) on the "momentum" 
plane: 

(23) n = {(m, n)} : 0<n<N, < m < M . 

The domain of -Fp,Q is the same, it is the Newton polygon for F(x 2 , y 2 ). On the boundary of the rectangular 
IT the eigenvalues of t m ,n as well as the values of Fp t Q are simple. Just putting e.g. y = in J^J, one gets 

(24) i ,oM) 2 -ti,o(x,0) 2 = (l-x 2 ) M . 

This equation defines all possible boundary eigenvalues t m fl. Subject of interest is the calculation of t m , n in 
the middle of II. 

Numerical calculations show that for all eigenstates the absolute values of t myn as well as the coefficients 
Fm,n grow significantly when (m, n) goes from the boundary of II to its middle. One eigenstate (up to the 
parity equivalence (|22|l ) is strictly separated from all others: absolute value of any its t m n is the maximal 
with respect to values of the same t m ^ n for all other eigenstates. We will call it the ground state. 

For a given eigenstate, especially for the ground state, the values of t m ^ n are maximal in some (m, n) — 
(Pq,Qo) in the middle of rectangular II. In the same point the coefficient Fp 0t Q has the maximal absolute 
value with respect to all other Fpq. The observed feature of the ground state is that Po+m '^ 0+ " with \m\ 
and \n\ being relatively small, depends essentially only on N/M and k when N and M are big. The same 

asymptotical independence of N, M is valid for ^ p "+ m 'Qo+" ag we \\^ Since Fp q takes the maximal value at 

f p ,Qo 



(Po,Qo)> expression l|17|) is the result of a competition between the domain of maximal values of Pp,q and 
big (or small) values of A MP /i N Q accompanying Pp,Q, si. (|13[1 . 

Another feature of {t m , n } may be mentioned. We observed that the sets of signs of {t m ^ n } for (m, n) 
surrounding (Po,Qo) are different (up to (|22|) ~) for different eigenstates. 

These observations allow us to suggest an idea for evaluation of l|15|l . Since both {i mj „} and {Pp,q} have 
a domain of dominance - the neighborhood of (Po, Qo) in the middle of II, far from the boundary - one can 
move to (Po, Qo) and concentrate on its neighborhood. The boundary of II is far from (Po, Qo), and in the 
limit N, M — > oo the boundary goes to infinity, so that the domain of the dominance becomes the open I 2 . 
It follows for finite TV, M, in the neighborhood of (Po + m,Qo + n) 

(25) \t Po + m , Qa +n\ 2 ~ \Fp a+m , Qo+n \ ~ e NM ^ K ^ , |m|and|n| are small , 

so that singular at N, M — > oo exponential is just the common factor for all eigenstates. Cancelling it, one 
does can evaluate the spectral equations l|15f) in the domain of dominance in the limit (J3J) . This will be done 
in the next section. 

5. Spectral equation in the thermodynamic limit 

N 

To rewrite equations or (|15|l in the thermodynamic limit N, M — > oo with £ = — being fixed, we 
need to introduce several notations. 
Define parameters c and a via 



(26) c = cot 



a jl + k 2 2 _ sin ^ 



2 V 3 - k 2 sin f 

At N, M — > oo the middle point (Po, Qo) is defined by 

(27) M ■ (l - = P - ui , N ■ (l - ~) = Qo - u 2 

where Po and Qo are even integers while u± and 112, —1 < ux^uz < 1, are fractional parts. If a is not a 
rational fraction of 7r, both u\ and 112 are extra variables. 

Define next the quadratic form parameterized in the terms of c and aspect ratio £: 

,™s ^/ \ 7T / 1 + C 2 9 1 - C 2 , 1 + C 2 9 

(28) n(p, q) = - [C^^P 2 + ^^P1 + C 

The AT, M — *■ 00 limit of l|15|) is based on the following behavior of the coefficients Pp,q: 

(29) Fp Q+p . Qa+q = (-)p+*+pi e NM *°^ ■ e-^P+u^+uz) . /q f l+ h + f2^{v + u 1 ,q + u 2 ) 



V 

where 

(30) So (k 2 ) 0(1, 1; k 2 ) = (l-g)log« 2 + 3<I>(|) -cj>^ 

Coefficients /o, /1, /2 are some functions of n 2 , (, u\ and 112 (a sketch derivation of l|29|) and the value of /o 
will be given in the appendix). 
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Define r TOi „ as the fine structure of i m ,„, 
(31) t Po+nhQa+n = yT^ hNMso{K2) r m ,n. 

Here, according to the idea of the previous section, we have moved to the middle (Pq,Qq) and canceled 
common exponential factor. Substituting (|29|) and l|31|) into (|15|) . cancelling the exponents and taking the 
limit N, M — *■ oo, we come to the following equations for r mj „, 

(oo\ \ ^ \ ^(\m+n+mn _ p -f2(p+ui ,q+u 2 ) 

\° ) / j / j \ I ' p+m.q+n ' p—m,q—n c 

The next substitution 

foo\ _ _ „-ifl(m+tii,n+ti 2 ) 

transforms 1|32|1 into the free from u\,U2 form: 

( 34 ) ^^(-) m+ " +mn e- n ("^c p _ m , 9 _„c p+m , g+ „ = 1 V P ,q e Z. 

m,nGZ 

Equations l|32[l and (|34[1 are two forms of in the thermodynamical limit. 

6. Analysis of 
For the analysis of l)34|l . let us modify it slightly at the first: 

where the cut-off parameter /3 > 1. 

Consider for a moment £ = 1. In this case 

(36) ft(m,n) = ^(c~ 1 {m + n) 2 + c(m - n) 2 ) , 
and we have two small parameters in (|35|1 . 

(37) g = e-' 37r/4c and Q = e"^ 4 . 

Equation (|35() may be analyzed in the terms of the perturbative expansion with respect to Q,Q. The zero 
order reads 

(38) <£, + o(l) = l => Cp, g = e M (1 + o(l)) , 
where e Ptq = (±) is the sign of c PtQ . In the first non-trivial order, 

(39) (l + ( 

This procedure may be continued, the result is a series with respect to Q and Q, 

(40) c p>g = e p , q (l + E X ( p m q n) Q m Q") ■ 



9 



with coefficients xp™ being sums of products of e m;n for (m,n) surrounding (p,q): The first few nonzero 
Xp™'"' with to + n < 4 are 

(41) Xp,g ^ = e P+l,9 e P— 1)9 "t" £ p,q— 1 £ P,9+1 ' 



(42) 

(2,2) _ , 
Xp.g — £p,g+l£p.g— l^p+l,g — l^p— l,g — 1 £p+l,g£p— l,q&p— l.g+l^p— l,q— 1 i ^p,g+l£p,g — — l,g+l 

~t~^p+l,g^p— l,g^p,g^p— 2,g ^p,g+l£p,g— l^p,g+2£p,g 1 ~t~ £p,q-\-\£p,q—\£p,q£p,q — 2 £p+l,g£p— l,g^p,g+l^p,g — 1 

l,g^p+2,g£p,g "I - £p+l,g£p— l,g^p+l,g+l^p+l,g— 1 7 

and 

(43) Xp 4 g 0) = ^p-l,g-l£p+l,g+l ) Xp ^ = £p-l,g+l£p+l,g-l ■ 

This procedure may be formulated for Q(m, n) with general £ as well. 

Conjecture 1. If [3 > 1, a/Z seria j^fj ) converge. Solution of is defined uniquely by the distribution of 
the signs e = {e m ,n}- 

Note, due to the parity structure of l|35f) . any distribution {e m . n } is equivalent to {e' m n = s m „(±) m (±)™}, 
sf. (EH). 

The homogeneous distribution e m ,n — (+) is the distinguished one since in this case c m ,n = cq, expression 
for co follows from (|35(l and matches the series form (|40|l . 

(44) Cq = I yj(_)ra+n+mn e -0n(ra,n) j 

\m,n / 

But if /3 — > 1, (@3J diverges as 

(45) c 



for some % = x( c ;0- This divergence may be explained by the term in 129JI : [3 = 1+ jy^.j w ben 

N, M — > oo , so that asymptotically 



'JVM 



(46) c = W-j— - y/NM. 

V /2X 

The distribution e m . n = (+) and c TOiTl = Co ~ V NM is the ground state according to the numerical tests. 
If the signs e m ,n vary for different to, n (even if only one sign is opposite to all the others), we have 

Conjecture 2. The seria \4-U\j with inhomogeneous e converge at f3 = 1. 
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We can explain Co ~ v NM in a bit different way. Consider for instance the following distribution of the 
signs: 



(47) 



(+) if n(m,n) < ~V, 



7T 

randomly (±) if fi(m, n) > —V 



In this case a very rough estimation gives 

(48) c Pi9 ~ W. 

Thus the finite- volume domain of positive signs on the infinite lattice is effectively equivalent to finite lattice, 
and the strip Q,(m, n) ~ — V plays the role of an effective boundary. 

The homogeneous distribution s m:n = (+) in a big volume V gives evidently the maximal eigenvalues of 
the quantum mechanical model, any variation of the signs gives an excitation of the spectrum. A distribution 
of the signs e mi>ni = £m 2 ,n 2 = ■■■£m h ,n k = (-) with (mi, rii)... (mfc, n*,) inside V and with all other e m , n = (+) 
inside V, is a candidate for a fc-particles state. 

One particle state, E mii „, = (— ) with all other £ m ,„ = (+) inside V, is described asymptotically by 
two continuous parameters = ( —=,—=). We expect a "dispersion relation" in the form 



<v vvj vv 

a smooth function of (/x, v) . The model evidently is gap-less. 

The behavior (|48|l allows one to suggest a candidate for the Hamiltonian of the system: 

(49) H = = -£<ne- n(m ' n) , 

sf. (|33|) . At the ground state w — /loV^, i.e. one can talk about the density energy — ho of the ground 
state, and the spectrum of H describes bound states — ho < — < 0. 

From the alternative point of view, one may consider the Hamiltonian 

(50) H' = —H . 

For this Hamiltonian, the ground state corresponds to a random distribution of the signs - we can say 
nothing about it. Excitations are the islands of constant signs in the sea of random ones, and its maximal 
value is described by the finite energy density +ho- 

7. Discussion 

The main results of this paper are the following. The ground state distribution of the moduli 

(51) r p . q = c p , 9 e-5 n (M) 
related to the moduli t m ^ n by i|31f) . is given by 

(52) r p , g = c e~i Q ^ , c ~ W 
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where V is the volume of the system. Any excited state is uniquely defined by a distribution of the signs 
£p,q = sign of Cp. q , the set of c p>q is the solution of 



The various problems of our approach are to be mentioned. At the first, we are still unable to evaluate (|53[) 
explicitly for non-periodical distribution of signs. Even f% in the asymptotic l)29|l or more exact estimation of 
H48fl are not known. Also it is not known yet how to express the momenta, corresponding to two orthogonal 
shifts of the periodical lattice, in the terms of {t m . n }. Without the physical momenta, one hardly can 
interpret physically the dispersion relation conjectured above. The third problem is the calculation of the 
spectrum of T{6\, 62, #3). It is still unclear how the transfer matrix of Zamolodchikov-Baxter model is related 
to {t m , n } for generic N, M . Some discussion repeating |13| is given in appendix. 
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Appendix A. sl 3 fusion algebra 

Let us demonstrate how the equations (|11I12|) generate the fusion algebra of auxiliary transfer matrices. 
Choose the particular value N = 3. The series Q may be rewritten as the four-term sum 
(54) 

N M 

J ^,y) = ^^(-i) m (-^T(-F y )\, n = 



n— m— 



(M \ / M \ / M \ / M \ 

£ {-xa x ) m t m , - (J2 0™ x ) m tm,i )yv y + ( £ (^rw )y 2 -(Yl H*O m *n»,3 y z ° v 
m=0 / \ro=Q / \ro=0 / \m=0 / 

= to(xo*) - h(xa x ) ya y + t 2 (xa x ) y 2 - t 3 (xa x ) y^o^ . 

One may show combinatorially jSj, tk{x) is the I4q=-i(slpf) transfer matrix for the Lax operators with the 
cyclic representation in the quantum space and the fundamental representation irk in the auxiliary space. It 
is supposed, ttq and ttn are the scalar representations, ~k\ is the vector representation etc. In N = 3 case 7T2 
is the co- vector representation. 

Decomposition of F(x 2 ,y 2 ) with respect to y 2 is following {to — e 27 "/ 3 and A 3 = a; 2 ): 

(55) F(x 2 ,y 2 ) = f[ ((1 - \^) M ~ V\l + ^ n ) M ) = A{x 2 ) - B(x 2 )y 2 + C(x 2 )y* - D(x 2 )y 6 , 

n=0 ^ ' 

where 

A(x 2 ) = (1 - x 2 ) M , D(x 2 ) = (1 + k 6 x 2 ) m , 

(56) ^--^{mr^-m M <-^ 



6„2\M 



C{x z ) = (1 + kV) 



1 — A \ / 1 — ALU \ I 1 — \Ll> 



1 + k 2 \J \1 + k 2 XloJ \1 + k 2 Xlv 2 
Equating now in all orders of y 2 , one comes at y° and y e to 

(57) t (x)t (-x) = (l-x 2 ) M , h{x)h{-x) = (1 + A 2 ) M 
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For the generalized chiral Potts model the choice is prescribed: 

(58) t (x) = (1 - x) M , t 3 (x) = (1 + \k 3 x) . 

The orders y 2 and y A give 

t^hi-x) = t a (-x)t 2 (x) +to(x)t 2 (-x) + B(x 2 ) , 

(59) 

h(x)t 2 {-x) = t 3 (~x)h(x) + t 3 (x)t 1 (-x) + C{x 2 ) . 
Relations (|59|l with (|58|l are exactly the fusion algebra for sis, H3 ED QUI ■ 

Appendix B. Asymptotic of Fp.q 

Let us discuss briefly the derivation of l|29|l . Taking into account l|2(JI) . one may use the saddle point 
method for the estimation of -Fp.Q. Basically, 

1 / / dXdY F(X, Y) 

(60) F ^ = w^ffic^^w- 

Let 

61 a„ = . a n = . 

Then 

(62) log (5^r) ~ NM (f( A > ^;^ 2 ) - v lo s A " v log " 

It has the extremum (minimum) with respect to A, \x (k 2 being fixed) at 1 

(63) r +r 2 = a p , r + r 3 = a q . 
The extremum value of g(A, fi; n 2 ) — — log A — — log /i is 

7T 7T 



(64) g(a p ,a q ;K 2 ) = - log* 2 + V $(r,) 

7T ^ — ' 

where the numbers 7"j are to be calculated via 

Oi + a 2 + a 3 a 2 + a 3 - a\ a 3 + a\ — a 2 ai + a 2 — a 3 

(65) r = 7T , n = , r 2 = , r 3 = 

and 

k 2 -1 

K 2 + l 



(66) a 2 = 7r — a p , a 3 = 7r — a g , ai = arccos I cos a 2 cos a 3 H — 2 - sin a 2 sin a 3 



sin t'q sin f\ 

The last equality is the solution of k 2 = with respect to a\. Therefore asymptotically 

sin r 2 sin r 3 

(67) Fp.Q = (- f+Q+ P Q . f Q . f\ + + ^ . pNMg( ap , aq ;^) 



^In details, A — = , fi 



ro + Ti as _ r + r 3 2 dg _ r 
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Function g(a q , a q ; k 2 ) has the maximum near a p = a q = ir — a, where a is defined by i|26|) . and 

1 + c 2 1 — c 2 
(68) g(a p ,a q ;n 2 ) = Q (k 2 ) - — (5a 2 + Sa 2 q ) - — Sa p Sa q , 



where Qq(k 2 ) is given by (|3U[) . Let further even integers Po, Qo and real numbers u\, u 2 are defined by l|27|) . 
Then 

(69) (5a p = — (p + ui), Sa q = —{q + u 2 ). 
Therefore, the leading term of 16711 is 

(70) Fp +P,Qo+ q = {-) p+q+pq ■ f ■ e NMao ^-^ p+u ^+ u ^ , 

where the quadratic form is given by (|28(l . 

The next order in (|67|l . F' = f\ + fz^Mj) + ui,q + U2), is the result of numerical tests. 

Appendix C. Theta-functions 

In the limit M,N — > 00 the polynomial F(X,Y) as well as the eigenstates of t a ,p(x,y) for periodical 
distribution of the signs s m , n become the theta-functions. In particular, equations 1)32(1 may be re-written 
in a theta-functions-like form: 

E2m 2n \ ( \ \n 2m+l 2n \ A \ ^ \m 2m 2n+l 

X V T 2m ,2n \ ~ I > ,(-) ^ 2/ T 2m+ l,2n I _ (Z^( _ ) 21 2/ ^2m,2n+l 



(71) 

" (H(~)" +m:E2m+1 y 2n+lr2 " l + 1 ' 2 »+ 1 ) = ^(-) p+9+M e-^ +Ul '« +U2 VV 9 
Let us re-define (x = e 11721 , y = e l1TZ2 ). Then the theta- function-like seria 



(72) T a<l3 (Zl,Z2) = (-) Qn+/3m ^n +Q ,2 n +/ 3 e l 



\-K{2m-\-Oi)z\-\-'\TT(2n-\- (3)nz2 

\ j i zm-\-a,~4n-\-p*~' 



stand for the transfer matrices. 

It is helpful to discuss some properties of theta-functions. Let 



(73) eW 2 (z 1)22 ) = ]Te 



P&1(p+ui ,q+U2 )+27rip2i +27riqS2 
p.q 



for our particular quadratic form £1 Q28JI. It has the general Jacobi transform property: 

(74) e^fe,^) = | e-^^—) egf^ (-«,,«!) . 

The other ^-function, related to F, is 

(75) F Uu u 2 (zi,Z2) = {-) p+q+pq e- n{p+Ul > q+U2)+2 ™ lP+2 ™ 2q . 
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One can easily see, 
(76) l 

F ui,« 2 ( z li z 2) — 2 (®«2«2 ( z l + 5; z 2 + ^) + 0«^,«2 ( z l + |) z 2) + 0li',«2 ( z li z 2 + ^) ~ ®ui,u 2 { z li z 2 )) 



For the case u\ = u 2 = the polynomial identity 

F 2 N,2M(x 2 ,y 2 ) = F N . M {x,y)F NM {-x,y)F NM (x,-y)F NM (-x,-y) provides 



4 ,z\ z 2sr _ ,Z\ + \ z 2 ^ i _ f Z\ z 2 + 1 N ,_ , #i + 1 + 1 , 



(77) /o/ r o,o( 2: i, ^2) — faFofl{ — 1 — )Fo,o( — ^ — iyKofY,- 1 ^ — )Fo,o(- , } 
The limit z\ , z 2 — > gives / for ||55J| : 



(78) /o = 

As well, the value of x for follows from 



f" 0,0(57 0)^0, o(0, 5)^0,0(2' 2) 



(79) ^ ( _ r+ „ +m „ e -^( m ,n) = r cj (0>0) = Ie£/«_eg> « (1-^ 

at /? -f 1 with x = eg + 2^M|^ =1 . 

Appendix D. Examples of periodical distribution 
Here we give an example is a periodical distribution of the signs. Let 

(80) £ 2m+ a,2n +/ 3 = £ Q ,/9 e^ um +^ 

with u,v — or 1. Periodicity of £ mj „ provides the periodicity of the series expansions (|40|) . and therefore 

(81) C2m+ ct) 2ri+/3 = £ 2 m+u, 2 n+f3 C-ufi ■ 

Equation gives 

for all four choices of (a, (3), its solution is Cq = c\ = Cq 1 = c\ 1 (it follows as well from the careful analysis 
of the structure of e^n-products in l|4l)|l). so that 



(83) c 2m+Q , 2n+/3 = e Q)/3 e i7r ( uro+w ' (e^-e^e^-e^e^-e^^eJI) 



-1/2 
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Appendix E. Transfer matrix of Zamolodchikov-Bazhanov-Baxter model 

In the last section we would like to describe the relation between (|20|) and Baxter's free energy for 
Zamolodchikov's model. We will refer to |12| . where the inhomogeneous model was considered and divisor 
parameterization was used. Equations (231) in l°°k like 

(84) J(X) T = T J(X') = . 

Here J{X) and J(X') are generating functions ©, operator T is a modified transfer matrix for 
Zamolodchikov-Bazhanov-Baxter model (in general, the Pauli matrices may be replaced by the Weyl al- 
gebra generators at root of unity). It follows from (|84|l . T up to a normalization is the product of algebraic 
supplements of J(X) and J(X'). 

In our particular case, J(X), J(X') and T after the quasi-diagonalization are 2x2 matrices (in the basis of 
the Pauli matrices). Transfer-matrix of Zamolodchikov's model T(0\, #2, #3), mentioned in the Introduction, 
is the trace of T: 

(85) T = Trace 2x2 T . 

Generating functions J(X) and J(X') stand for J(X(X) N / 2 , fJ,(X) M / 2 ; k 2 ) and J{X{X') N / 2 , f i(X') M / 2 ; k 2 ) in 
the present notations, where 

/ SR x 2 . 2 i shifts sin ft 

(86) K = tan — = - — - — : — — 

2 sm ft sin ft 

is the ^-parameter in both J(X) and J(X'), and explicit evaluations for A and /i from |12j to the terms of 
linear excesses ft give 

(87) X(X) = e-^+W^ , rtX) = e'^' + ^^- , 

sin P3 sin P3 

and 

(88) X(X') = e'^'+M^- , n{X') = e-'^+M^ . 

sin P2 sin P2 

It gives us the identification {rj} = {a permutation of f3j} and relates (|20|l to Baxter's answer for the 
partition function per site k: 

(89) logfc = normalization + ^ (|l log 2 sin ft +$(/?,) ) . 

The reader may see the discrepancy, —i- log 2 sin /3j in (|89|l and — log 2 sin ft in Q2()jl. it means that the 
normalization is not trivial - it comes from a certain variational principle. 
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